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We consider a small multi-polaron model obtained by coupling the many-body Schrodinger equation for N 
interacting electrons with the energy functional of a mean-field crystal with a localized defect, obtaining a 
highly non linear many-body problem. The physical picture is that the electrons constitue a charge defect 
in an otherwise perfect periodic crystal. A remarkable feature of such a system is the possibility to form 
a bound state of electrons via their interaction with the polarizable background. We first prove that a 
single polaron always binds, i.e. the energy functional has a minimizer for A'^ = 1. Then we discuss the case 
of multi-polarons containing N > 2 electrons. We show that their existence is guaranteed when certain 
quantized binding inequalities of HVZ type are satisfied. 
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1. Introduction 

A polaron is a quantum electron in a polar crystal which is able to form a bound state by using 
the deformation of the medium which is generated by its own charge [Tj. Likewise, a multi-polaron 
or N -polaron is the system formed by the interaction of N electrons with a polar crystal. 
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That a polaron can be in a bound state is a rather simple physical mechanism. When the 
(negatively charged) electron is added to the medium, it locally repels (respectively attracts) the 
other electrons (respectively the positively charged nuclei) of the crystal. A local deformation is 
thus generated in the crystal, and it is itself felt by the added particle. In other words the additional 
electron carries a "polarization cloud" with it. It is therefore often useful to think of the polaron 
as a dressed particle, that is a single (composite) particle with new physical properties: effective 
mass, effective charge, etc. For an A^-polaron the situation is a bit more involved. Since the effective 
polarization has to overcome the natural Coulomb repulsion between the particles, bound states 
do not always exist. 

The question of what model to use to describe the polaron is an important and non trivial one. 
In the Born-Oppenheimer approximation, a quantum crystal is a very complicated object, made of 
infinitely many nuclei and delocalised electrons. The accurate description of such a system is a very 
delicate issue and, for this reason, simple effective models are often considered. They should remain 
mathematically tractable while still capturing as much of the physics of the system as possible. 

A famous example is the model of Frohlich [8i dating back from 1937, in which the crystal 
is described as an homogeneous quantized polarization field with which the electrons interact. In 
the limit of strong coupling between the electrons and the field, the model reduces to Pekar's 
theory [211 [22j [23l [171 [20] . There the crystal is a classical continuous polarizable model, leading to 
an effective attractive Coulomb interaction in the energy functional of the theory: 

^ JR3 ^ JR3 Jk3 \x - y\ 

Here tp is the wave-function of the electron, em is the static dieletric constant of the crystal and 
we work in atomic units. The variational equation corresponding to (|l.ip is sometimes called the 
Schrodinger-Newton or Choquard equation. 

It is the attractive Coulomb term in (jl.ip that leads to the existence of bound states of elec- 
trons, i.e. minimizers (or ground states) of the energy functional. Whereas the energy functional 
for electrons in vacuum has no minimizer, Lieb [15^ proved the existence and uniqueness (up to 
translations) of a ground state for Pekar's functional (|l.ip . 

The same nonlinear attractive term is obtained in Pekar's model for the A'^-polaron. Then, as we 
have already mentioned, depending on the strength of the attractive Coulomb term as compared 
to the natural repulsion between the electrons, one can get binding or not. It is an important issue 
to determine in which paramater range binding occurs |10[ [HI [3 [T3] . 

The approximations made in the construction of Frohlich's and Pekar's models reduce their 
applicability to situations where the iV-polaron is spread over a region of space much larger than 
the characteristic size of the underlying crystal. One then speaks of large polarons. In |14j we have 
introduced a new polaron model by coupling the energy functional for electrons in the vacuum to 
a microscopic model of quantum crystals with defects introduced in [21 [3] . Unlike in Frohlich and 
Pekar theories we take the crystal explicitly into account and make no assumption on the size of 
the electron. Our approach thus qualifies for the description of both small and large polarons. The 
model takes the following form (for one electron): 




Here V^^^ is the (periodic) electric potential generated by the unperturbed crystal, which is felt by 
any particle added to the system. The nonlinear effective energy i^crys represents the interaction 
energy between the electrons and the crystal. It is defined using a Hartree-Fock theory for the 
response of the electrons of the crystal to a charge defect. The state of the Fermi sea of the 
perturbed crystal is given by a one-body density matrix 7, that is a non-negative self-adjoint 
operator on L'^{R^). As in [2l [3] , we write 

7 = 7per + Q (1-3) 
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where 7pp;. is the density matrix of the periodic unperturbed crystal and Q is the local deformation 
induced by the charge defect \ip\'^. The effective energy F^ys then takes the form 

Fery.[|VP]= inf (f [ ^4*^dxdy + J-„y.[Q]y (1.4) 

Three main ingredients enter in 1)1.4^ : 

• Electrons are fermions and must thus satisfy the Pauli exclusion principle, which gives in 
the formalism of density matrices the constraint < 7 < 1 as operators. This justifies the 
constraint on admissible perturbations Q imposed in (|1.4p . 

• The electrons forming the polaron interact with the perturbation they induce in the Fermi 
sea. This is taken into account by the first term in p.4p where pg is the charge density 
associated with Q, given formally by pq{x) = Q{x,x) (we use the same notation for the 
operator Q and its kernel). 

• Generating a deformation of the Fermi sea has an energetic cost, represented by the func- 
tional J-"crys in (|1.4|) . The somewhat complicated definition of this functional will be recalled 
below. It was derived in [2]. 

More details on how we arrived at the form above can be found in the introduction of |14j . Let us 
mention that this model only takes into account the displacement of the electrons of the crystal 
and neglects that of the nuclei. This is arguably an important restriction, but our model already 
captures important physical properties of the polaron, and on the other hand this is all we can 
treat from a mathematical point of view at present. 

In this paper we will show that a (single) polaron described by the energy functional ()1.4|) always 
binds. The case of iV-polarons is more sophisticated, as now the effective attraction resulting from 
the polarization of the crystal has to overcome the electronic repulsion. The energy functional 
corresponding to ()1.2|) in the case of the iV-polaron is given by 
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\'^{xi,...,xn)\' , 

-off[*J = / I n 7 J V:r, X-ATjl i- ^ j j \dXl---dXN 

\Xk — Xf,\ 



l/p",,p* + F„.y,[p*] (1.5) 



where p>p is the usual density of charge associated with the many-body wavefunction ^ whose 
definition is recalled in ()2.16|) below. 

In fact, our model (|1.2p is closely related to Pekar's functional. We proved in [14] that Pekar's 
theory can be recovered from (|1.2|) in a macroscopic limit where the characteristic size of the 
underlying crystal goes to 0. In addition to clarifying the physics entering the Pekar model, this 
result also gives some interesting insight on the model (|1.2p . in particular, regarding the question 
of the existence of binding. Indeed, it is known [T3] that Pekar's functional has a ground state in 
some range of parameters. We deduce in [14] that sequences of approximate minimizers for (II. 5p 
converge in the macroscopic limit to a ground state of the Pekar functional, thus showing that our 
model at least accounts for the binding of large polarons in this regime. In this paper, we want to 
derive conditions ensuring that there is binding in the case of small polarons where the macroscopic 
limit argument and the link to Pekar's theory are irrelevant. 

Quite generally, for many-body quantum systems, the existence of bound states of N particles 
depends on the validity of so-called binding inequalities. If E{N) denotes the infimum energy of 
some physical system containing N particles, a ground state containing N particles exists when 

E{N) < ^mm^E{N ~ k) + E°^{k) (1.6) 
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where E°°{N) denotes the energy of the same N particle system, but with all particles 'sent to 
infinity'. For example, for atoms or molecules comprising N electrons, E{N) includes the contri- 
bution of the electric potential generated by the fixed nuclei, while E°°{N) does not. Particles 'at 
infinity' no longer see the attraction of the nuclei. Note the formal similarity between inequalities 
(|1.6p and those appearing in Lions' concentration compactness principle |18| I19j . an important 
mathematical tool used in nonlinear analysis. The major difference is that the former are quan- 
tized and thus more difficult to relate to one another. See (13; for a more precise discussion of this 
connection. 

It is not difficult to discuss on physical grounds why inequalities (|1.6|) are sufficient for the 
existence of bound states. Indeed, (|1.6|) says that sending particles to infinity is not favorable 
from an energetic point of view. In mathematical terms, the inequalities ()1.6p avoid the lack of 
compactness at infinity of minimizing sequences. The existence of a ground state then follows from 
the local compactness of the model under consideration. Nevertheless, the mathematical proof that 
inequalities of the type (|1.6p are sufficient for the existence of bound states of A^-particles is highly 
non-trivial because the problems E{N), E{N — k) and E°°{k) are set in different Hilbert spaces. 
In the case of atoms and molecules, the fact that inequalities of the form (|1.6p imply the existence 
of bound states is the content of the famous HVZ theorem, first proved independently in [TH \^ 

In this paper we prove an HVZ-type theorem for our polaron functional (jl.5|) when N > 2. We 
have to face two difficulties. First the functional is invariant under the action of arbitrarily large 
translations (those leaving invariant the periodic lattice of the crystal), so the energy functional 
does not change when particles are sent to infinity. The correct binding inequalities therefore take 
the form 



E(N) < min E(N - k) + E(k). (1.7) 

k=l...N-l 

Second, the energy contains the highly nonlinear term Fcrys[p*]- We are thus faced with the com- 
bination of the difficulties associated with many-body theory and those inherent to nonlinear 
problems. A general technique has been introduced in |13) to tackle these questions. Our purpose 
in this paper is to explain how one can deal with the model ()1.5|) using the method of 13 . Our 
main task will be to control the behavior of the (highly nonlinear) effective polarization energy 

^ crys ■ 
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2. Statement of the main results 
2.1. The mean-field crystal 

We begin by recalling the precise definition of the crystal functional entering in ()1.2|) . More details 
can be found in [2 El E] . 

We fix an ^-periodic density of charge ^p^^ for the classical nuclei of the crystal, with ^ 
a discrete subgroup of R'^. It is enough for our purpose to assume that ^p^^ is a locally-finite 
non- negative measure, such that /p /ip^,,. = Z E N, where F — M.'^/^ is the unit cell. 

In reduced Hartree-Fock theory, the state of the electrons in the crystal is described by a one- 
particle density matrix, which is a self-adjoint operator 7 : L^(M^) — L^(R'^) such that < 7 < 1 
(in the sense of operators). When no external field is applied to the system, the electrons arrange 



Binding of small polarons 



5 



in a periodic configuration 7 = 7pci., which is a sohition of the reduced Hartree-Fock equation^ 

'7^er = l(-oo,ep)(-A/2 + Vpy, 

- AVpO,, = 47r(p^o_^_ - , (2.1) 



Here pA denotes the density of the operator A which is formally given by pa{x) = A(x^x) when 
A is locally trace-class. Also, l(„oo,eF) ('^) denotes the spectral projector of H onto the interval 
(— oOjEf)- The real number ep in (|2.ip is called the Fermi level. It is also the Lagrange multiplier 
used to impose the constraint that the system must be locally neutral (third equation in (|2.ip ). The 
unique solution to the self-consistent equation ()2.ip is found by minimizing the so-called reduced 
Hartree-Fock energy functional [SJ [5] • 

We are working in the so-called atomic units in which the mass m and the charge e of the 
electrons of the crystal are set to to = e = 1. Also we neglect their spin for simplicity (reinserting 
the spin in our model is straightforward). 

By Bloch-Floquet theory [21], the spectrum of the ^-periodic Schrodinger operator 

is composed of bands. When there is a gap between the Zth and the {Z -\- l)st bands, the crystal 
is an insulator and ep can be any arbitrary number in the gap. As in [5], in the whole paper we 
will assume that the host crystal is an insulator. 

Asumption 2.1 (The host crystal is an insulator). 

The periodic Schrodinger operator H^^^ has a gap between its Zth and {Z + l)st bands, and we fix 
any chemical potential ep in the corresponding gap. 

When the quantum crystal is submitted to an external field, the Fermi sea polarizes. The 
method used in [2] to define the energetic cost of such a polarization relies on the following idea. 
The energetic cost to move the electrons from 'y^^^ to 7 is defined as the (formal) difference between 
the (infinite) reduced Hartree-Fock energies of 7 and of 7pcr- Denoting by 



Dif, g) := I I MMdxdy [ dk (2.2) 



the Coulomb interaction (where / denotes the Fourier transform of /), one arrives at the functional 

T,,ys[Q] Tro {{H^,, - £p)Q) + ^DipQ^pQ) (2.3) 

where Tro denotes a generalized trace, see (|2.6p and p.9p below. For convenience we also denote 

Tcrys[p,Q] -.^Tio {{Hi, -eF)Q) + \D{pQ,pQ)+D{p,pQ). (2.4) 

The functional setting in which the terms of these equations make sense is defined as follows. Any 
operator Q satisfying the constraint 

- 7pcr < Q < 1 - 7per (2-5) 

is decomposed as 

Q = Q— + Q-+ + Q++ + Q+- (2.6) 



^Sometimes called Hartree equations in the physics literature. 
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where Q = 7porQ7pcr; Q ^ = TpcrQ (l ~ 7por) ; ^^^^ SO On. It is proved in ^ that for Q satisfying 
(|2.5p and v e L^(K'^) n i^(K'^), FcTys[v, Q] is finite if and only if Q is in the function space 

Q = {q e &^\q = Q\ |V|Q e Q^+^Q^' e 6\ |V|g++|V|, |V|Q— |V| e e^} (2.7) 
that we equip with its natural norm 

IIQIIq = IIQIIe^ + ilQ++llei + \\Q''\\& + ll|V|g||e^ + II |V|g++|V|||ei + |||V|Q— |V|||ei. (2.8) 

The symbols 6^ and & denote the Schatten classes of trace-class and Hilbert-Schmidt operators 
on L^(]R^) respectively (see [IS] and [H]). For operators in Q, the kinetic energy in (12.31) is defined 
as 

Tro(ff°,, - £f)Q = Tr (|i/p"„, - eFr/'(Q++ - Q-)\HI,, - eF\^'^) , (2.9) 

see [2]. More generally, one can define the generalized trace as 

TroQ = TrQ+++TrQ"- (2.10) 

when and Q^^ are trace-class. Note that Trg differs from the usual trace Tr, the operators 
in Q not being trace-class in general. They nevertheless have an unambiguously defined density 
Pq E LI^^{R'^) (see [2], Proposition 1). It belongs to L^(M^) and to the Coulomb space 

C = {p\D{p,py^^ (2.11) 

and it holds by definition 

Tro(Fg) = / VpQ (2.12) 



for any V eC. 

Having defined in (|2.3p the total energetic cost to go from 7p(,j. to 7p(,j. + Q, we can give a sense 
to the energetic response of the crystal to an external density i^. The state of the Fermi sea is 
obtained by solving the following minimization problem 



FcrysW]= jj^i^ „ {Div,PQ)+Tcrys[Q] 



(2.13) 



As shown in [2], for any i/ g L^(R'^) n _L^(R'^), this minimization problem has at least one solution 
in Q. The corresponding density pq is in L'^{M.^) but in general it has long range oscillations which 
are not integrable at infinity [1]. 

2.2. The small polaron 

To our crystal we now add N quantum particles, which are by assumption distinguishable from 
those of the crystal. In reality they are electrons having the same mass m = 1 as those of the 
crystal, but we want to keep m arbitrary to emphasize that in our model the additional particles 
behave differently from those of the crystal. This will also allow us to compare with the results we 
have obtained in [14]. 

The total energy of the system now includes the term -FcryslH with v = \'ip\'^ (polaron) or = p^, 
(iV-polaron). For the single polaron, the energy is given by 

I f:5^|VV'(:c)|'+C(^)l^(^)l') ^^ + ^".ys[|^H. (2.14) 
For the A^-polaron with > 2 it reads 



^{xi,...,xnW 

l<k<e<N 

V;°„(x)p*(a:) dx + Fcrysip*]. (2.15) 
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As we think that there is no possible confusion, we do not emphasize the particle number N in our 
notation of the energy £. The density pq, is defined as 

p<ii{x) — N j \^(x,X2, ■ ■ ■ .xn)]^ dx2 ■ ■ -dxN- (2-16) 

JR3(N-1) 

The corresponding ground state energies read 

E{1) := inf |f [^], ^ e i/i(K3), |^|2 = i| (2.17) 

and 

£;(7V) := inf * e i?i(M3"), ^' fermionic, / l^-p^ll. (2.18) 

Here by 'fermionic' we mean antisymmetric under particle exchange: 

*(xi,...,Xj,...,Xj,...,a;Ar) = -^[xi,...,Xj,...,Xi,...,XN) for any i (2-19) 

as is appropriate for electrons. Recall that we have neglected the spin for simplicity. 

We now state our main results. In the single polaron case we are able to show the existence of 
a bound state. 

Theorem 2.1 (Existence of small polarons). 

For N = I, we have 

E{1) < i?pe. := inf a (-^A + V^^,,^ . (2.20) 

There are always minimizers for E(l) and all the minimizing sequences converge to a minimizer 
for E{1) strongly in H'^{W^), up to extraction of a subsequence and up to translations. 

Inequality (|2.20p expresses the fact that binding is energetically favorable : the right-hand side 
is the energy an electron would have in absence of binding. 

In the iV-polaron case we can give necessary and sufficient conditions for the compactness of 
minimizing sequences. 

Theorem 2.2 (HVZ for small iV-polarons). 

For N >2, the following assertions are equivalent: 

(1) One has 

E{N) < E{N -k) + E{k) for allk^l,...,N-l. (2.21) 

(2) Up to translation and extraction of a subsequence, all the minimizing sequences for E(N) 
converge to a minimizer for E(N) strongly in H^{M.^^). 

Remark 2.1. For this result, the fermionic nature of the electrons is not essential. The same 
theorem holds if the electrons are replaced by bosons , i.e. the wave function 5" is supposed to be 
symmetric under particle exchange. 

As discussed in the introduction, this theorem is rather natural from a physical point of view. It 
is not expected that the conditions (|2.2ip hold in general. As in Pekar's theory, one should expect 
the existence of minimizers to depend on the choice of parameters entering the functional (in our 
case only the periodic distribution /itpg^ of the nuclei). Testing the validity of these inequalities is 
a challenging task that would require more knowledge on the properties of the crystal model than 
we presently have. In particular, the decay at infinity of the minimizers of the crystal model should 
be investigated. 

In [14j we have considered a macroscopic limit in which V^^^ is replaced by m~^V^^^{x / m) and 
-Fcrys[p*] is replaced by m~'^ Fc-cys[m? p^^im-)]. By scaling, this is actually the same as letting the 
mass m of the polarons tend to zero in the model considered in this article. In the limit to — >■ 
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the ground state energy Em{N) converges to Pekar's energy involving the macroscopic dielectric 
constant Em of the crystal defined in [4] (up to a simple oscillatory factor, see [l^ for details). It 
was shown in [13j that the binding inequalities are satisfied in Pekar's theory when em is large 
enough. We conclude that in this case they will also be satisfied for m small enough and therefore 
minimizers do exist in this case. 

The rest of the paper is devoted to the proof of Theorem 12.21 One of us has considered in 
Section 5 of j l3| a general class of nonlinear many-body problems of the form 

/ \\y^\^x,^{xi,...,XN)\^ + V \^{xi,...,XN)?W{xk-xi)\dxi---dxN^F[p^] 

and provided sufficient assumptions on the interaction potential W and the non linearity F under 
which a HVZ type result similar to Theorem l2.2l holds. The assumptions on W include the Coulomb 
interaction we are concerned with in this paper but, unfortunately, our crystal functional Fc^ys 
does not seem to satisfy all the properties imposed on F in [13] . Also the presence of the periodic 
potential V^^^ adds a new difficulty. Nevertheless the general strategy of |T3] still applies and our 
goal in this paper is to explain how to overcome the difficulties associated with -Fcrys- 

Section [3] gathers some important properties of the crystal functional that are to be used in the 
proofs of Theorems 12.11 and 12.21 presented in Sections H] and [5] respectively. 

3. Properties of the crystal energy 

In this section we roughly speaking prove that i^crys satifies Assumptions (Al) to (A5) of |13j . 
Section 5. We are only able to prove a little less, but the properties we do prove are sufficient for 
the proof of Theorem 12.21 as we explain in Section [5l 

We start in Section 13.11 with almost immediate consequences of the defintion of -Fcrys, and 
devote Section 13.31 to the more involved fact that our crystal functional satifies a 'decoupling at 
infinity' property. The proof of this property requires some facts about localization operators that 
we gather in Section [321 



3.1. Concavity, subcriticality and translation invariance 

The following is the equivalent of Assumptions (A4) and (A5) in [TF, Section 5. 
Lemma 3.1 (Concavity). 

-Fcrys is concave on {p Cz C, p > 0}. Moreover it is strictly concave at the origin: 

F„y,[ip] > tF,,ys[p] (3.1) 

for all p Cz C, p >0 and all < t < 1. 

Proof. The functional J-"crys [p, Q] defined in (j2.4p is linear in p. As by definition 

= inf {Tcrysip, Q], -7pcr < Q < 1 " 7pcr} , 

it is clearly a concave functional of p. As for the strict concavity we note that 

J"crys[ip, Q] = Tro - ep) g) + \d{pq, pg) + tD{p, pq) > tT,,y,[p, Q] > tF,,y,[p] 

for all < i < 1 by positivity of the kinetic and Coulomb energies. Taking for Q the minimizer 
corresponding to tp which is known to exist by [2] l4l proves p.l|) . □ 

The next lemma will be useful to prove that minimizing sequences for our polaron model are 
bounded in H^(R.^^). It is the equivalent of Assumption (A3) in fT^, Section 5. 

Lemma 3.2 (Subcriticality). 

The functional Fays is locally uniformly continuous on L^^^ . More precisely, we have 

|F„ys[p] - Fcrys[p']| <C\\p- p'Wler. (3-2) 
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for a universal constant C > 0. Moreover, for every e > 0, we have 

> F„M'] >-e [ |V<^|2 --( f 1^ 

for all ipe H^{R^). 



, 3 
|2 ' 



(3.3) 



Proof. For any p € L^/^ and any Q € Q we can complete the square in the electrostatic terms of 
•^crys [p, Q] and obtain 

J-„ysb, Q] = Tio - ep) Q) + ^DipQ +p,pQ+p)- ^Dip, p) > -^Dip, p). 

Taking the infimum with respect to Q and applying this with p = \(p\'^ immediately yields 

FcrysM] > -lD{\ip\\ > -C Mi,.,, 

by the Hardy-Littlewood-Sobolev inequality ([E], Theorem 4.3). Using now the Sobolev and Holder 
inequalities we get as stated 



Then, replacing p hy p — p' we also have 
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J-e,y.[p - p', Q] > --Dip -p',p- p'). 



Choosing now for Q a minimizer of J^^rys [p, Q] we deduce 

i^crys[p] - F„ys[p'] > -^D{p- p',p~ p'). 

Without loss of generality we can assume that the left-hand side is negative. We conclude that 
there exists a constant such that 

\Fcrys[p\ - Fcrys[p'] \ < C \\p - ^11^6/5 

using the Hardy-Littlewood-Sobolev inequality again. □ 

Finally, we note that our functional is invariant under the action of the translations of the 
periodic lattice .if. Note that in [13], full translation invariance is assumed (see Assumption (A2)). 
However, what is really used in the proof of the results there is the invariance under the action of 
arbitrarily large translations. 

Lemma 3.3 (Translation invariance). 

For any p G L^/^ and any translation t £ ^ of the periodic lattice, 

F,,y,[pi-+T)]^F,,y,[p]. (3.4) 

Proof. We denote by Q a minimizer of Fcrys[p, Q]- Clearly pq{- + f) = putQUf where Uf is the 
unitary translation operator acting on _L^(R'^) and defined by Uff = /(• — r). We deduce 

FcrysM- + r)] < J-„.ys[p(- + f),^QT] = Tro (r (i/p",, - ep) f*Q) + \d{p, pq) - D{p, pg) = i^„ys[p] 

by translation invariance of the Coulomb interaction and the fact that -ffp^r commutes with the 
translations of the lattice Exchanging the roles of p(--f r) and p and applying the same argument 
proves that there must be equality. □ 
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3.2. Some localization properties 

In order to prove that the crystal energy of two distant clusters of mass decouples we will use a 
localization procedure. We here provide several facts about this procedure that will be useful in 
the next section. 

We will use the following local compactness criterion in Schatten classes. Its standard proof is 
omitted. 

Lemma 3.4 (Local compactness in Schatten spaces). 

Denote by & the class of compact operators A of some Hilbert space f) such that {TT{\A\P)f^P < 
+00, with the convention that &°° denotes the class of compact operators. 

• If An A weakly-* in and K, K' e 6°° then KA„K' KAK' strongly in 6^ 

• If An A weakly in 6"^, K <E & and K' £ & then KAnK' KAK' strongly in & with 
l/s = l/p+l/q + l/r. 

We next prove an interesting property concerning the localization of functions of the Coulomb 
space C. 

Lemma 3.5 (Localization in the Coulomb space C). 

Let X € L°°{MP) be such that Vx G £^(R'^). Then XP G C for all p E C. More precisely, there exists 
a universal constant C > such that 

IIxpIIc<^^(IIxIIl» + I|VxIIl3)IHIc- (3-5) 

Proof. Recall that C is the dual of the homogeneous Sobolev space II^{M.^). For any V G C^(IR."^) 
we have by the Sobolev inequality 

XPV <(47r)-i/2||p||J|V(xV^)||^. 

< (47r)-i/2 WpWc ( WxWl^ l|vy||^. + llVxILs \\v\\^, ) 

<(47r)-i/^ \\p\\c{\\x\\L^+C\\Vx\\L^)m\\^. 
and the result follows. □ 

Note that for a typical smooth localization function defined by xr{^) — xi^/ R) with x = 1 
in -8(0, R), X = in -6(0, 2), the right-hand side of (13.51) does not depend on R: 

||V(xKf)HL / xl\'^XBf = 2'R-'f x{x/Rf\Vx{x/R)\'dx^2^ [ x'|Vx|' (3.6) 

This is what makes the lemma particularly useful. For such a localization function we also have 
the following approximation result: 

Lemma 3.6 (Approximation in the Coulomb space). 

For any p E C, we have {xrYp ^ P strongly in C. 

Proof. Recall that ||XflP||c — ^IIpIc "^ith a constant C that is independent of R, by (j3.6p and 
Lemma rS. 51 So we can use an "e/2 argument", and replace p by a function p' e L^/^, for which 
the statement follows from the Hardy-Littlewood-Sobolev inequality. □ 

It will be important in the sequel to know that weak convergence of a sequence (Qn) in Q 
implies strong local compactness of the corresponding pq^ : 

Lemma 3.7 (Local compactness for charge densities). 

Let (Qn) be a bounded sequence in Q such that Qn — ^ Q weakly in Q. Then pq^ pq weakly in 
L^ C\C and strongly in L\^^. 
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Proof. We know from Proposition 1 in 2 that pQ^ pg weakly in L^nC. Only the strong local 
convergence is new. The latter follows from the statement that 

xQnX -> xQx 

strongly in the trace-class for every smooth function x of compact support. To see this, we 
write as usual 

Qn^Q++ + Q-+ + Qt- +Q-- (3.7) 

and consider only the first two terms as the other two can be dealt with in a similar way. We have 

xQt^X = {x (-A + 1)-^/^ }{ (-A + 1)^/^ (-A + 1)^/^ }{ (-A + 1)-^/^ x}. 

The operator x(^A + 1)^^/^ is compact and (—A + A + 1)^/^ converges towards 

(—A + l)^/^(5++(— A + 1)^/^ weakly-* in by assumption. By Lemma 13.41 we deduce that 
xQt^X XQ^^X strongly in 

We argue similarly for the off diagonal terms, writing this time 

xQt-x - {x (-A + 1)-^/^ }{ (-A + 1)^/^ Qt-}{i",..x}- 

Again the operator x(^A + l)^^/^ is compact and we can write 

7perX = 7per(ffpcr + M) (^^per + A^)"' (1 - A) (1 - A)-^. 

Here ^ is a large enough constant such that Hp^,. > ~p/2. The operator '-fp^j.{Hp^j. + p) is bounded 
by the functional calculus. Also, (i?por + m)^^ (1 ^ A) is bounded by Lemma 1 in '2|. Finally, 
(1 - A)-^x G 6^- Thus 7°^^% G Since (-A + 1)1/2q+- ^ (_a + 1)1/2q+- weakly in 6^ by 
assumption, we deduce by Lemma [3.41 again, that xQi x ^ xQ^ X strongly in (3^. 

We have proved that xQnX ^ xQx strongly in 0^, but PxQ^x ~ X^Pq„i so we deduce that 
PQ„ PQ strongly in L\^^. U 

Finally, the following algebraic property, which does not seem to have been noticed before, will 
be very useful when constructing trial states for the crystal functional. 

Lemma 3.8 (Adding states using localization). 

Let n he an orthogonal projector on a Hilhert space S}, and Xj V self-adjoint operators on 
such that -l- < 1 • We introduce the corresponding localization operators 

a: = n^n + (1 - n)x(i - n) and r==n?7n + (i-n)7;(i-n). 

Let Q, Q' two self-adjoint operators such that — H < Q, Q' < 1 — H. Then we have 

- n < XQX + YQ'Y < 1 - n (3.8) 

as well. 

Proof. We have ~XUX - YHY < XQX YQ'Y < X{1 - n)A + Y{1 - I1)Y . Now, the result 
follows from the estimate 

xnx + YUY = nxHxn + nr^nr^n < n(x^ -f ?7^)n < n 

and the equivalent one 11 replaced by 1 — 11. □ 

Typically, one can think of x and r] as being two localization functions satisfying x^ + = 1 ■ 
In our context, the projector H will be j^^^. and the main virtue of the corresponding operators X 
and Y is that they commute with 7°^^. 
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3.3. Decoupling at infinity 

Here we provide the most crucial ingredient of the proof of Theorem 12.21 namely the fact that the 
crystal energy of the sum of two distant pieces of mass is almost the sum of the energies of these 
pieces. This is the content of the following proposition, which is the equivalent of assumption (A3) 
in [13]. Note however that we prove much less than what is stated there. Fortunately, the proof of 
Theorem 25 in [T3] does not actually require such a strong assumption as (A3), as we will show in 
Section [5] below. 

Proposition 3.1 (Decoupling at infinity). 

Let {pn) be a bounded sequence in the Coulomb space C such that pn p weakly. Then 



In the above, one should think of p„ as being constituted of two clusters of mass, p and p„ — p, 
whose "supports" are infinitely far away in the limit n — > oo. This is mathematically materialized 
by the weak convergence to of p„ — p. The proposition then says that the total energy is the sum 
of the energy of the pieces, up to a small error. Proving p.9p is a difficult task because of the long 
range behavior of the response of the crystal : it is known [3] that the polarization pq of the Fermi 
sea has long range oscillations that are not integrable at infinity. The oscillations generated by p 
are seen by p„ — p (and conversely) but, fortunately, they contribute a small amount to the total 
energy, which is controlled by the Coulomb norm and not the norm. 

Assumption (A3) in [13j is a little different from (13. 9|) . There it was assumed that Pn — Pn + Pn 
where and are bounded in L^/^ and that the distance between their supports goes to infinity, 
with no assumption on the size of these supports. In Proposition 13. II it is implicit that one of the 
two clusters of mass has a support of bounded size and is approximated by its weak limit p. This 
additional assumption is harmless for our purpose because we are dealing with a locally compact 
problem. 

In the course of the proof of Proposition 13. II we will establish the following, which we believe 
is of independent interest. It gives the weak continuity of the (multi-valued) map p„ i— >■ Q„ = 
argmin TaysiPn, ■]■ 

Corollary 3.1 (A weak continuity result for J^ays)- 

Let (pn) be a bounded sequence in the Coulomb space C such that pn p weakly and Qn be any 
minimizer of TcrysiPm ■]■ Then, up to extraction of a subsequence, Qn ^ Q weakly in Q where Q 
minimizes Tcrys[p, ■]■ 

We now present the 

Proof of Proposition [37ll We begin with the more difficult part, that is the proof of the lower 
bound corresponding to p.9p : 

Step 1: Lower bound. We denote by Qn a minimizer for Q t-^ J-crys[pn, Q]- Corollary 2 in [2 states 
that the energy functional J^ciys[PmQ] controls the norm ||Q||q- 



The upper bound is obtained by taking a trial state Q = 0. Using that p„ is bounded in C, hence 
that D{pn, Pn) is bounded, we deduce that the sequence (Qn) is bounded in Q. Up to extraction of 
a subsequence, we can assume that Qn ^ Q and, by Lemma 13.71 that pq^ pg weakly in nC 
and strongly in Ll^^. 

We now introduce a smooth partition of unity + 77^ = 1 such that x = 1 on the ball -6(0, 1) 
and X = outside of the ball B{0, 2). Similarly, = 1 on \ b{0, 2) and 77 ee on B{0, 1). We 
also assume that Vx and V?7 are bounded functions. Then we introduce xr{^) ■= x{^/R) and 




(3.9) 



> J^ciys[pn,Qn] > C\\Qn\\Q - -D{p„,Pn)- 
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rjii{x) = r/{x/R). In our proof R is fixed and will go to infinity only in the end, after we have taken 
the limit rt — > oo. Since 

is independent of R and ||X-rIIloo ^ 1 ^-iid ||?yi?||^oo < 1 by definition, we conclude that 

||pq.'7h|Ic + II^Q^XflJIc < C' (3.10) 

uniformly in n and i?, by Lemma [3.51 Also, {xb.)^Pq„ ~^ {xrYpq strongly in and weakly in 
L^, thus also strongly in L^/^ by interpolation. In particular, {xrY PQ^ ~^ {xr)^PQ strongly in C 
by the Hardy-Littlewood-Sobolev inequality. 
We start by writing 

D{PQ^,Pn) = D{PQ„,P) + D{{VrT PQn., Pn - P) + D{{xRf PQ,„ Pn - p) ■ 

Note first that 

lim D{pQ^ , p) = D{pQ, p) and lim D{{xr)'^Pq„ , p„ - p) = 

since /o„ — p ^ weakly in C and {xrYpQtl ~^ {xr)'^PQ strongly in C. So we conclude that 

lim (d{pq^, p„) - D{pQ, p) ~ D{{r]R)^pQ^, p„ - p)) =0. 

Then we write in a similar fashion 

D{PQ^ , PQ„ ) = DiixRfpQ^ , {xr)^PQ,. ) + 2i?((xfl) Vq„ , iVR^PQ^ ) + Vq„ , (riRfPQ. ) 

and conclude that 

Jim (z?(pQ„,PQj-i?((xfl)VQ, (xi^.)VQ)"2D((xi^)VQ, (^fl)VQ)-^((^7fl,)VQ„, ('7i?)VQj) = 0. 

Following ideas from [TTl [5] , we now consider a better localization method for Qn which respects 
the constraint (|2.5p . We define the two localization operators 

^R = 7porXi?.7por + (l - 7pcr) XR (l ^ 7pcr) 

Yr = 7pcr^m7pcr + (1 - 7pcr) VR (l " 7pc,) (3.11) 

that have the virtue of commuting with the spectral projectors -fp^j. and (1 — 7pgj.). Note that in [2], 
the choice Xr ~ -^/l — is made. Here we change a bit the strategy and we only have 

Using Lemma 10 in [2], it can actually be shown that + > 1 — C/R for R large enough. 

As noticed first in |llj , the main advantage of these localization operators is that they preserve 
the constraint. Simply, using that Xr-i^^^Xji = 7perXfl7pcrXfl7pcr < 7pcr(Xfl)^7pcr < 7per and the 
similar estimate Xr{1 — ^p^^)Xr_ < 1 — 7por: see that when — 7per < Q < 1 ^ 7pcrj then 

-7per < ~Xr-1%,Xr < XrQXr < Xr{1 ~ J°p,,)XR < (1 - 7°^). 

These localization operators will be used to localize the kinetic energy, i.e. write it as the sum of 
the kinetic energies of XrQuXr and YRQnYR plus errors. As we have already used the localization 
functions XR and t^r to deal with the Coulomb terms, we also have to relate PXrQXr to PxrQxr 
and PyrQYr to PriRQriR - These two tasks require the following lemma, whose proof will be detailed 
in the Appendix (it is based on several estimates from [2]). 

Lemma 3.9 (Properties of the localization operators Xr and Yr). 

There exists a universal constant C > such that 

Tro(ir°er " £f)Q - Tro(i?;^,, - eY)XRQXR - TvoiH^.r ~ eF)YRQYR 



<15^\\Q\\q, (3.12) 
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WxrQxb^ - XrQXrW^ + WmQm. - YnmRlQ < ^ IIQII 



(3.13) 



and 



WPx 



c 



for all QeQ. 

Using p.l4p and the facts that Qn is bounded in Q and /9„ — p is bounded in C, we get 

D{Pr,RQ^VR^Pn - P) - D{pYrQ„Yr, Pn " P) 

where C is independent of n. Similarly, using p. 141) and (I3.10[) . we deduce that 

T^nR ; PvrQtiVr ) I D{pYnQ^Yii , PYrQ„Y„ 



(3.14) 



C 
< — 
- R 



C 
< — 
- R 

Now if we use (I3.12p to deal with the kinetic energy, we arrive at 

-^crys[Pn] ^ CYys{Pn:Qn] 

= Tro(iJ^er - eF)XRQnXR + Tro(iJ°er - £F)YRQnYR + D{pQ,p) + D{pYrQ,^Yr. Pn - p) 

+ \D{{xRfpQ, {XR?Pq) + D{{xRfpQ, imfpg) + ^D{pYRQ„YR,PYRQ^Yn) + £n{R) 

> Tto{H^,, - eF)XRQnXR + F,,y,[pn -p]+ D{pQ,p) 

+ \D{{xRfpQ, {xnfpQ) + D{{XR?PQ, {mfpQ) + £n{R) 

where we have used that YRQnYR is an admissible trial state for J-"crys[pn — and where 

C 



limsup \en{R)\ < 



(3.15) 

with a constant C that is independent of n. Passing to the liminf and using Fatou's lemma yields 
liminf (Fc,y,[p„] - F„y,[p„ - p]) > Tyo{H^o, - £f)XrQXr + D{pQ,p) 

+ ^DiixRfpQ, {xr?Pq) + D{{xR?PQ. {mfpQ) - ^ 



Using Lemma 13.61 we have 

lim D{{xRfpQ,{xRfpQ) =D{pq,pq) 

R^oo 



and 



lim D{{xRyPQ.{mYPQ)=^- 



Using Fatou's lemma again for the kinetic energy term, we arrive at the result 

liminf {F„ys[Pn] - i^cry«b„ " p]) > "TioiH^,, - £f)Q + D{pQ, p) + ]-D{pq,pq) > (3.16) 
which is the lower bound corresponding to 



Step 2 : proof of Corollarv \3.1\ with p = 0. 

We pick a sequence p„ ^ and denote by (Qn) the corresponding sequence of minimizers. 
Since (pn) is bounded in C, (Qn) is bounded in Q and, up to extraction, converges weakly to some 
Q ^ Q, which implies that pQ^ pq weakly in C. We prove here that (3 = 0. 

Thanks to the lower bound part of p.9p we have just proved, we write 

-l^D{pQ, pQ)+F„ys[pn]+o{\) < Fc,ys[-PQ]+Fcrys[Pn]+o{l) < i^crys [Pn-Pg] < FcrysW]- D {PQ , PqJ , 

using Qn as a trial state for Fciys[pn — Pq] and the simple lower bound J\;rys[i^, Q] ^ ~'^D{i',i'). 
Taking the limit n — > oo we therefore obtain D{pq,pq) — 0, which implies Q = 0. 
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Step 3: Upper bound. We now construct a trial state for Fcrys[pn] to obtain the upper bound part 
of p.9p . In previous works [H] [2], the special structure of the set of admissible states was used 
(see the Appendix of [11]). We propose here a new method based on Lemma [3^ 

Let Q and Qn be two minimizers for, respectively, the problems i^crys[p] and -Fcrys [Pn — p] ■ Recall 
that they must satisfy the constraint — 7per ^ Qj Qn < l~7per ^nd that, using Step 2, Qn 0. Let 
Xb. be a localization function of compact support as before, G N, f G and rj^ = xr{- ^ Nf). 
Consider the trial state (we use the notation of Lemma \3M with H = 7p(,r) 

where C/jv? is the translation unitary operator, which commutes with H^^^. hence with jp^^. When 
N is larger than twice the diameter of the support of XRi then xji + Xr{' ~ -^^) ^ 1 ^^d therefore 
we have — 7p(,r — Qn.r. < 1 ~ 7por' by Lemma We deduce, using Lemma 

^l:rys[p«] < ^crys[Pn, Qn,b] = TrQ ( (7per " *^f) XrQXp;) + ^(pn, Pq) + ^D{xrPq, XrPq) 
+ Tro ((7pcr - £f) YjiQnYfi) + D{pn - P,PQ„) + ^D{T]lpQ^,T]jiPQj 

+ D{xrPq,VrPQ,J + D{p,i]jiPQj - D{7]%pQ,pn) - D{pn - P,XrPqJ + £n{R) 

< i^crysW + i^crys[p« - P] + D{XrPQ,VrPqJ + D{p, r]%PQ,J - D{l]jiPQ,Pn) 

- D{pn - p, XrPqJ + £n{R) 

where e„(i?) satisfies (13.15^ and we have used the fact that x'rtVr: ^r: < 1 and the translation 
invariance of the energy. Now, X?jPQ„ ~^ strongly because of the strong local compactness of pg^ 
(^Lemma 13.71) . Using also p„ p, pq^ and Lemma l3.6l we can take first the limit n — > oo and 
then the limit i? — > cx) to conclude 

hmSUp {Fcrys[p„] - Fcrys[p] ~ Fcrys[Pn " p] ) < 
n— ^oo 

and the proof of Proposition 13.11 is complete. 

Step 4: End of the proof of Corollary \3.1\ Let {pn) be any sequence such that pn ^ p weakly in 
C, and Qn be any associated sequence of minimizers for J-"crys[Pra, •]• Extracting a subsequence we 
may assume that Qn ^ Q in Q. Coming back to the lower bound p.l6p obtained in Step 1 and 
using (|3.9p . we see that 

-?^crys[p] = lim [FcrysiPn] " -^crys [Pri - P] ) > -^crys[P, Q] > -FLrys [p] • 
n— >oo V / 

This shows that Q minimizes J^:;rys[p, •] and concludes the proof of Corollarv l3.1l □ 
4. Existence of small polarons: Proof of Theorem 12.11 

Before turning to the more complicated case of N particles for which we have to adapt Theorem 
25 in |13] , we deal with the simpler one-particle case. The proof that a minimizer always exists for 
one particle follows from usual techniques of nonlinear analysis. In this context the most difficult 
is to verify the one-particle binding inequality (|2.20p . which we do first. 

Step 1. Proof oj the one-particle binding inequality 

The aim of this first step is to prove the following important 

Lemma 4.1 (One-particle binding). 

We have 

Eil) < Ep,,. inf a A + y^o \ (4^^^ 
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Proof. Let Wpor denote the first ^-periodic eigenfunction of —A/ {2m) + V^^^., which is a solution 
of 



A , 



(4.2) 



We assume that Wpor is normahzed, /p luporP — 1 where F is the unit cell of Since Uper € 



-ffpgj.(r), we have also Vpc 



^pcr| 



G iJp(,j.(r). The Fourier coefficients {i^pcr {k))ke^* thus satisfy 



(4.3) 



(|fcp9per(fc))fee^* G and consequently belong to ^^(if*) : 

|t'pcr(fc)| < oo. 

fee.Sf* 

Here ^* is the dual lattice of whose unit cell will be denoted by F*. We can write 

l"per(a^)l' = ^ I] Ppor(A:)e''=-" 

Consider now a fixed function x G C^(R"^) such that J |xp = 1, and define the following test 
function for the variational problem E{1): 



ipx ■■= Upcr{x)x\{x), with x\ix) := A ^^^x (^) • 
The corresponding density is 

|Va(x)P :=|upe.(a:)P \xxix)\' = ^ ^ i>per(fc)e^'- \xxix)f 

Remark that 



(4.4) 



(4.5) 



D{\xx\'e^''-,\xx\'e^'-)=in\- 



■ dp 



47r 



|p/A+fc|2 ^ A^oo A3|fc|2 

for any fc ^ 0. Using (|4.3p and the fact that Upcr is normalized, we deduce that 

/ 1 



x\Hp) 



dp 



\\i^x?-\xx\\ = o 



VA3/2 



Similarly, the normalization factor is 



\Uperix)\'^\Xxix)\'^ dx 



(27r)3/2 



for all p e N. Of course, we have by scaling 



D{\xx\',\xxn = jD{\x\'Axn 



We deduce from all this that 



crys 



/m3 Hxl 



= Fcrys[\Xxn+0^-^^ 



by (|3.2p . In Theorem 1.4 of [M] we have studied in detail the behavior of the crystal energy when 
the external density is very spread out. We have proved that 



F^rys [IXAp] = F,,y, [A-3|x(7A)n = [Ixl^ 



(4.6) 



where is Pekar's effective interaction energy 
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Since Em > 1, we have [p] < for all p. So the exact (first order) behavior of the crystal energy 



for our trial state is 



^ crys 



/r3 IV'. 



A ' V A 



The two other terms in the energy £ are easier to handle. A simple computation based on the 
equation (j4.2p of Upor shows that 



Jr3 2m ^ 7r3 2m Jr3 

(see Lemma 2.2 in [14]). Of course, 

since Wpcr € ^^por C L°°(R'^). As a conclusion we have shown that 




^■A ,^SlM1,„ 

— -C/por H 1- Oa^-oo I 



'/r3 IV-A 

Since [\x?] < 0, the inequality (|i?T|) follows. □ 
S'fep ^. Compactness of minimizing sequences and existence of a minimizer for N — 1 



We now turn to the proof of the other statements in Theorem 12.21 dealing with the one-particle 
case £'(1). 

Let (ipn) be a minimizing sequence for £'(1). Since > £crys[|'0P] > —D{\ip\'^, |i/'p)/2, it is easy 
to see that {ipn) is bounded in H^{M.^). We define the largest mass that subsequences can have up 
to translations by 

M := sup I / IV'P : 3(a;fc) C R^ V„,(- - Xk) ^ weakly in H\R^) 



We know [18] that M = if and only if tpn strongly in LP{R^) for all 2 < p < 6, a 

lie 



phenomenon that is usually called vanishing. But if this is the case, we get l^nnL — > by 



the Hardy-Littlewood-Sobolev inequality, and therefore £crys[|V'n| ] — >■ by p.2|) . We then get 
E{1) > Eper '■= inf f (-f^pcr) which is impossible by Lemma H?T] Thus Af > 0. 

Since Af > we can find a subsequence (denoted the same for simplicity), such that tpni'^Xn) 
•0^0. We can of course write Xn = kn + Un where fc„ £ ^ and y„ € F. Extracting subsequences 
again we get Un ^ y & ^- Therefore tpni' — kn) tpi- + y) ^0. Since our energy functional is 
invariant under the translations of the new sequence ■(/'«(• — kn) is again a minimizing sequence. 
Without loss of generality we can thus assume that ipn '4' 0- Now, if we can prove that 
/r3 IV'P = 1) ■'^6 '^ill g6t strong convergence in and it is then standard to upgrade this to strong 
convergence in H^. We argue by contradiction and assume that < /jj3 I^Z-P < 1. 

We will now show that the energy decouples in two pieces. Since ipn ^ 4' ^'^ H^{M.^) we may 
assume that IV'nP ^ IV'P in C We then use that, by (|3.9p in Proposition 13. 1[ 

£cry«[|V'«|'] > F„ysM^] + F,rysUn\^ " +o(l). 

Note that 

- I^P - - V-P = 23?^(V„ - V) ^ 

strongly in L-'^(]R'^) (we use here that ipn strongly in Lf^^ and an e/2 argument), hence in 

^^'^ by interpolation. Thus 



lim 

n— ^oo 
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by ()3.2p in Lemma 13.21 and we arrive at 



^crys 

On the other hand, it is clear from the weak convergence ^pn ~^ V' iii (and from the fact 

that the form domain of -A/(2to) + V°^^ is H^{R^}}, that 



Hence we have shown that 



£[i^n] > £[i^] + £[i^n - i^] + 0{1). 



Now we use that Fd-ys is concave to infer 



crys 



/k3 \tpn - V'l^ 



leading to 



/r3 IV'n - V'l^ 



o(l) 



> fiV-] + ^JV« - ^Pj £^(1) + 0(1). 
Passing to the limit n — > oo, we find 

i?(i)>f[^]+(^i-^jV|2)i?(i). 

It is now time to use the strict concavity at the origin p.ip 

Fcrys > (^^ I^P) i^c,y> ^ 

which yields 



/r3 IV'I^ 



/k3 1^1' 



> 



iV'ni^(i) 



Therefore we have proved that £'(1) > E{1) which is a contradiction, unless J^^ \ip\'^ = 1. This 
concludes the proof in the case of one particle. 



5. Binding of small A/'-polarons: Proof of Theorem 12.21 

We now turn to the case of > 2. With the imput of Section |3l the proof more or less follows that 
of Theorem 25 in [13]. We nevertheless sketch the main steps for the convenience of the reader. 
We will denote 

j=i ^ ' i<j 

In order to relate problems with different particle numbers to one another, it is crucial to introduce 
the antisymmetric truncated Fock space 

N n 

J-<A'^0/yL2(R3) 

n=0 i=l 
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where /\ is the antisymmetric tensor product and we use the convention L^(]R'^) = C. A state 
on J"-^ is an operator T G 6^{J^-^) with Tr(r) = 1. In the sequel we restrict ourselves to states 
commuting with the number operator 

AT 
n=0 

This means (see [l^, Remark 7) that they take the form 

T = Goo® ■■■OGnn (5.1) 
with Gu e el (ALi i^(K^)). We denote by 

N 

H = 0i7(n) 

n=0 

the many-body second-quantized Hamiltonian. To any state F are associated a density p-p S ^^(R'^), 
one-body density matrix [rj^'^ e &^{L'^{R^)) and two-body density matrix [r]^'^ g G^[L'^(]^3^ x 
L^(R^)) (see [T3], Section 1). We can extend the energy to Fock space as 

£[r] ^Tr^<«(Hr)+Fc,ys[pr] 

= Tri2(K3) ((-A + V^pO,,) [r]i-i) + Tri.(R3)>,i.(R3) (^[r]^^^) + i^crys[pr] 

where W acts on L^(]R^) x L^(R^) as the multiplication by \x — y\^^- For a pure state F = 
® . . . ® I*) (^"1 with * e L^(R^^) one can check that £[F] More generally, for a state of 

the form ()5.ip . we have 

AT 

£:[F] = TiA? L^(R-^) (ff (")G„„) -I- F, 

n=l 

Step 1. Large binding inequality. 
We claim that 

E{N) < E{N -k) + E(k) for ah fc = 1, . . . , - 1. (5.2) 

To see this, we consider the following trial state: 

■^^ := ^^-'^ A (. - Rf) (5.3) 

where (^^~^) and (vE''^) are compactly supported fixed trial states for E{N — k) and E{k) respec- 
tively, T G .if is a lattice translation and i? € N is large enough for p^N-k and p^t (. — Rt) to have 
disjoint supports. The symbol A denotes the antisymmetric tensor product. We first take the limit 
i? ^ cx) to obtain 

E{N) < £[^^-''] (5.4) 

Optimizing then with respect to vj/^-'^ and ^I''^' concludes the proof of ()5.2p . To see that (|5.4p holds, 
we note that by construction 

Pij,N = pq,N-k + p-q,k (. — Rt) 

for large enough R, thus we can use Proposition 13.11 and take the limit R oo with R E N to 
obtain 

lim Fciys[p^N] = -Fcrys[p*w-fc] + -Fciysb*'-']. 
it— ^oo 

The other terms in the energy can be treated as usual to obtain ()5.4p . 

Note that the argument here also proves by contradiction that Item (2) of Theorem 12 . 2 1 implies 
Item (1). If there is equality in (|5.2p . we can choose VE'^"'^ and 'i'^ minimizing sequences for 
E{N — k) and E{k) respectively and, taking i?„ — > oo very fast, we obtain a minimizing sequence 
for E{N) that is not precompact, even up to translations because some mass is lost at infinity. 



N 
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Step 2. Absence of vanishing. 

We consider a minimizing sequence (5'n) for E{N) and denote by r„ = ® ® |5'„) the 

associated state in the truncated antisymmetric Fock space. It is easy to see, using in particular 
Lemma [X^ that (5'„) is bounded in H^{M.^^). As in the one-body case treated before, we define a 
criterion for the vanishing of the minimizing sequence. We use the concept of geometric convergence 
(see Section 2 in [T3] for the definition). We fook at the the mass of the possible geometric limits, 
up to translations and extraction, of (r„) 

M := sup {Tr (MT) , 34 C R^ 4r„X- T} 

where we recall that Af is the number operator in Fock space. As explained in [13], Lemma 24, 
if M = then p^^ strongly in L^'(R'^) for all 1 < p < 3. Using then Lemma [3.21 we obtain 
-Fcrys[/3*„] and therefore 

E{N) = lim Sl^n] > Ma(H{N)] = NEp^, 

71— >00 \ / 

where 

^(^) - E {-^ + - E (K^X. • (5-5) 

Note that, by induction on A^, (|5.2I) implies E{N) < NE{1). We have already seen in (j4.ip above 
that E{1) < inf a (iJpor)- Hence we reach a contradiction and conclude that M > 0. 



Step 3. Decoupling via localization 

Since M > (and arguing as in the previous section) we have, up to the extraction of a subsequence, 
Vri^nV^ L with Tr(A/'r) > and where (z/„) C .if is a sequence of lattice translations. Using the 
invariance of the energy. Lemma 13.31 can thus assume that our minimizing sequence satisfies 

r„ r (5.6) 

with Tr(A/'r) > 0. Also we have y^prv ^ yfpr weakly in 7J^(R'^) and strongly in Ly^^. Also 
pr^ Pv in the Coulomb space C and we immediately deduce by (j3.9p that 

-F'crys[/3r„] > -Fcrys[/3r] + -F'crysbr,, - Pt] + o(l). 

We now pick a sequence of radii i?„ — oo and define smooth localization functions X-R„ ^-nd 
r]R^ such that Xr^ + Vr„ = 1, supp(xfl„) C S(0,2i?„) and supp(x/?„) d M? \ B(0,3i?„). For 
any bounded operator B (in particular the multiplication by a function x) on L^(M'^) such that 
< BB* < 1 we will denote by (r)^ the B-localization of a state F, as defined in [TSl, Section 3. 
Of importance to us will be the following properties of localization: 

Pr^ = X^Pr 

[T^Y-' = xn-\ 

[r^]''' = x®x[r]'''x®x- (5.7) 

Also, for a state of the form (j5.1l) . writing 

(r)x/?.„ — Gq"" ffl • ■ • ffl G^''" , (F),,^^ = Gg''" ® . . . ® G^" , 
the condition x\„ + = 1 implies the relation 

Tr(Gf")=Tr(G:^"^.). (5.8) 

Using concentration functions as in Step 4 of the proof of ^12]) Theorem 25 we have, extracting 
a further subsequence if necessary 

(r»)xH„ ^ r strongly in (J"^^) (5.9) 
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and 

(Xfl„)Vr„ ^ Pr strongly in LP{R^) for all 2 < p < 3. (5.10) 
Using p.2p . this can be used to prove that 

^crys [pr„ - Pr] = F„ys [(?7i?„)^/0r„] + o(l). 

Thus 

^crys [Pr„] > Fays [Pf] + ^crys P(r„),^ + o(l). 

We have seen that the nonlinear energy fcrys decouples. The other terms are treated follow- 
ing [13]. For the one-particle part we use the IMS formula 

to obtain (we use (15.71) ) 



2m ' 'p-y'-"^^"J ; i?2 

= Tr (H(r„)^,J + Tr (H(r„)„,J - — . 

The Coulomb interaction is treated exactly as in |13| and we conclude 

(^'„,i/(iV)^'„) > Tr(H(r„)x„ J +Tr(H(r„)^„J +o(l). 
Using Fatou's lemma as well as the strong convergence of (Xfl„)^Pr„, we finally get 

>£[T]+£ [(r„)^Hj + o(l). (5.11) 
which is the desired decoupling of the energy. 

Step 4- Conclusion. 

The rest of the argument follows exactly [13] . Writing the geometric limit of r„ 

F = Goo ffi . . . © Gat TV, 



and using the concavity of i^crys, the fundamental relation (|5.8I) as well as the convergence (|5.9p . 
we arrive at 

N 

E{N) > ^ Tr(G,,) {E{j) + E{N - j)) . 

Assuming the strict binding inequalities (j2.21l) . this is possible only when Gn = . . . = Gn-in~i = 
0. Hence we necessarily have Gnn 0, otherwise we would obtain a contradiction with the fact 
that Tr(A/'F) > 0. 

To conclude, it is then enough to prove that Goo = Oi which is an easy consequence of the strict 
concavity p.ip of F^rys (see Step 5 of the proof of Theorem 25 in [TB] for details). We deduce that 
Tr(GArjv) = 1 Tr(|*„) (*„|), hence that the weak-* convergence of |*„) (*„| in 6^ {L'^ {R^))to 
Gnn is actually strong because no mass is lost in the weak limit. As Gnn — \'^) (^| where ^' is 
the weak limit of 'i/n, we conclude that converges to ^I' strongly in L'^{M.'^). The convergence in 
iJ^(M^) follows by standard arguments. □ 
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Appendix. Proof of Lemma 13.91 

We follow ideas of [5]. We start with the second estimate p.l3p of Lemma 
instance 

{XRQXR-XRQXRy 



and write for 



= XiiQ--Xr ~ XrQ'~XR - bpcr; Xfl]QXK7por " Xfl7porQ [XR. 7por] 

= {Xr - xr)Q''Xr + XrQ''{Xr - xr) 

- [7pcr> XR]QXRlpcr - XRlpcrQ [XR ]. (5.12) 

When we multiply by |Hp(,,. — eFl""^^^ on both sides, we have to estimate several terms. For the last 
one we write 



per 



1-1/2^0 



epl --7pcrl^^pcr - ^f\'/'Q [Xfl„7pcr]|ffpcr " ^^1'/' 



+ XR^U^per - EFr/'Q [XR, IpeMpo 

7° 111 
/pcrJ II q2 



1/2 



(5.13) 



/perl per ■~r\ LA.-fi.i /pcrJI per 

In Lemma 10 of [2] it is proved that || [xr^ 7per] ||p;2 = 0{R~^) and the same proof can be employed 
to show that 



[XB.i Tper] l-^per 



,1/2 



as well. Also, following the proof of Lemma 11 of |2l, one can show that 



eF^/^Xi^l|i^pcr-eF| 



-1/2 



For the first term of the right-hand side of (I5.13P we can thus write 



1^, 



per 



I per I per 



pcrJ I per 



< 



l-ffpor - £f| 



1/2 



XR 



\Hpc 



£f 



-1/2 



j'\H'-e^\'/'Q ^„ [xR..^tML--eF\'/' 



< CR- 



We turn to the second term. Using that \H^f,^ — £f|^/^(1 — A) and its inverse are bounded by 
Lemma 1 in [2], we see that 

|i/"-£F|l/2Q| <C 



e2 



(l-A)i/2g <C||g||g. 



So we get 

\xRlpjKor - £f|'/'Q [Xfl,7pcr]l^^p°cr " ^F 

Then, for the first terms in (j5.12p . we use that 
|H°,, - SFl'/'iXn - xr)Q-Xr\H%, - e^\''^ 

< C WQWq |||7?pV - ep\'/^iXn - xr)\H"p,, ~ e^r'/^ 
Of course we have 



1/2 



< 



SI 



R 



|-H"per ~ Sf\^^'^Xr\HI^^ - £f| 



-1/2 



l-ffpcr - <^f\^^'^Xji\H° - Ef 



-1/2 _ 



\HIc^ " ev\^^'^ ,Xb^ |iJp„. -£f| ^^"^ + Xf 



per ^^1 i^^it l-'-^pcr 

and a similar formula for Xr ~ xr- The first commutator is a 0{R^^) by Lemma 11 in ^2j. Also 
we have \\Xr — xr\\ = 0{R^^) by Lemma 10 in [2]. In total, we have actually shown that 

\Hl, - e^\'/'{XRQXR ~ xrQxrVIK.^. - £Fr/^ < 
In order to conclude the proof of p.l3p we still have to show that 



R 





per 



EFr^iXRQXR-XRQXR] 



S2 



< 



g||QII( 

R 



(5.14) 
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The proof of this is similar but simpler and we omit the details. 

By Proposition 1 of [2], we know that the linear map Q E Q t-^ pg E C O is continuous. So 
the third estimate (|3.14p in the statement simply follows from the estimate p.l3p on the norm of 
Q. 

It remains to discuss the kinetic energy estimate p.l2|) in Lemma 13.91 First we remark that 



(Xfl.) = IpcrXR (7pcr + (7per)^) Xfl7pcr = (^fl.) + [7pcr, Xb] {lpor)^[XB. l\ 



I per 

and a similar equality for {x%)^^ • Since by construction 
this yields 

ixl)- + (Xfl)++ - ^fl, - [7per, XB? and 
From this we deduce that 

Tro(i?°,, - £F)g = Tr(i/0,, - ep){Q++ + Q-) 

= Tr 



Tpcr] 



{vl)- + iv'B)^^=Y^-h%VR] 



= Tr 



ixl + VBXKor - ^f) + (gp°e.. - SfKxI + rf^) 



(Q++ + Q— ) 



[7p°or, XB]HKor ' ^f) + (gp"e.. - £F)[7ge., XR? 
_ hper. m]\H^e. ' ^f) + (gger ' £F)[7gor: ^7^1' 

hence that 

Tro(ff° , - ev)Q - Tro(i?° - eF)XRQXR - TvM^, ~ ev)YRQYR 



1 



Tr {[Xr, [Xr, HlJ] + [Yr, [Yr, HlJ]) {Q++ + Q 



Tr ([7°cr, '7fl]'(^^p°cr - Ef) + (i^pV. ~ eF)[7pcr, m?) + Q"")- 



We conclude that 

|Tro(iJ°3, - eF)Q - Tro(ii;^,, - £F)Xi?Q^fl, - Tro(i7;^,, - eY)YRQYR\ 

<C\\Q\\a ( \Hl, - e^\-^'^[XR,[XR,HlJ]\Hl, - e^\-^'^ 



+ 





per 



,|l/2 



The last two error terms can be estimated by arguing as before, using (|5.14p . For the double 
commutators, a computation shows that 

[Xr, [Xr.HIJ] = (7;^cr)^([Xfl,7peJ [Xfl, A] + [Xfl, A] [Xfl,7pcr] + |VXflp)(7°or)^ 

- 7per([Xfl,7pcr] [Xfl., A] + [Xfl,, A] [Xfl,7pcr] + I VXiJp) 7pcr- 

We have [xk, A] = (Axk) + 2iVx • p with p = — iV. Using then that p|^^per ~ EfI^^^^ is bounded 
and the fact that || [Xfl; 7pcr] || = 0{R~^), we conclude similarly as before that 

\HpcY - £f\^^^^[Xr, [Xr,H°]]\H°. - EfT^^^ 
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The term involving is treated similarly. This ends our proof of the first estimate p.l2|) , hence 
the proof of Lemma 13.91 □ 
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